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Abstract. We consider by means of Monte Carlo simulations the relaxation in the 
paramagnetic phase of the anti-ferromagnetic Ising model on a triangular lattice and 
of a fully-frustrated Ising model on a square lattice. In contradistinction to previous 
studies of the second model, we show that spin-spin correlation functions do not decay 
with a stretched-exponential law at low temperature but that both models display an 
exponential decay with logarithmic corrections that are interpreted as the signature of 
topological defects. 
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1. Introduction 

The study of relaxation in frustrated systems is of particular interest because their dy- 
namics may become anomalously slow at low temperature and eventually freeze below 
a certain temperature. The paradigmatic example of such a slow dynamics is given by 
spin glasses which combine both frustration and randomness. Frustration alone does 
not necessarily imply a glassy dynamics. Even though each one of its plaquettes is 
frustrated, the anti-ferromagnetic Ising model on a triangular lattice [2] (AFIM) was 
shown to display the same dynamics as unfrustrated systems [17]. Interestingly, if elas- 
tic deformations of the lattice are allowed and coupled to the spin degrees of freedom, 
the relaxation is well described by a (stretched-exponential) Kohlrausch- Williams- Watts 
law like in glasses [IJ[6j[T5]. This behavior is observed experimentally for example with 
closely-packed colloidal spheres [7J. 

In this paper, we consider both the AFIM and the Fully-Frustrated Ising Model 
(FFIM) defined on a square lattice by the Hamiltonian 

- f3H = ^ [(Tx,y<7x+l,y + (-l/^^j/^jz+l] , Vx,y = ±1 (1) 

x,y 

where f(x,y) = x + y in the so-called zigzag model and f(x,y) = x in pileup-domino 
configuration. In both cases, each plaquette of the square lattice contains an odd number 
of anti-ferromagnetic bonds and, as a consequence, is frustrated. The pileup-domino 
configuration allows for exact diagonalization of the transfer matrix by free fermion 
techniques [12J. Both AFIM and FFIM are known to belong to the same universality 
class. At the critical temperature T c = 0, spin-spin autocorrelation functions decay 
algebraically 

C(t jS ) = (o-(t)o-(s))~(t-sy^ (2) 

where 77 = 1/2 and z = 2 is the dynamical exponent. In the case of the AFIM, evidences 
have been given of the existence of topological defects interacting through a logarithmic 
Coulombian potential in the paramagnetic phase [HI EB] • While no signature of these 
defects was observed in earlier simulations of the AFIM [9], we have shown that such 
defects manifest themselves in the FFIM [13]. When the system is initially prepared in 
the paramagnetic phase and then quenched at T c = 0, they pin the domain walls and slow 
down their motion so that aging takes place in the same way as in homogeneous systems 
but with logarithmic corrections. The correlation length grows as £(£) ~ (i/lnt) 1 / 2 and 
spin-spin autocorrelation functions behave as 

where both z and A are compatible with 2. 

This result suggests that at finite temperature, the relaxation should be exponential, 
like in unfrustrated systems, but with logarithmic corrections due to the presence of 
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topological defects. However, on the basis of Monte Carlo simulations, a stretched- 
exponential relaxation [j] 

C{t) ~ e- (i/r)K (4) 

with a temperature-dependent exponent k has been reported in the FFIM below the 
temperature T p ~ 1.701 at which Kasteleyn-Fortuin clusters start to percolate [2j|3]. In 
this work, we reconsider the AFIM and the FFIM and show that Monte Carlo data are 
in better agreement with an exponential decay when taking into account logarithmic 
corrections than with a stretched exponential. In the first section, a dynamical scaling 
hypothesis is set up to predict the behavior of equilibrium spin-spin autocorrelation 
functions in fully-frustrated Ising models. The expression is then compared with Monte 
Carlo data for the AFIM in the second section and the FFIM in the third one. 



2. Dynamical scaling of spin-spin autocorrelation functions 

We analyze the two-time spin-spin correlation functions in the framework of dynamical 
scaling [Hj. Upon a dilatation with a scale factor b, the equilibrium correlation 
C(r, t, T) = (<Xo(0)ay(i)) at temperature T is assumed to satisfy the homogeneity relation 

C(f, 1/t, \T -T c \) = b~ 2x °C(r/b, b z /t, \T - T c \b x ' v ) (5) 

where x a is the scaling dimension of magnetization density with 2x a = 77 for two- 
dimensional systems and z is the dynamical exponent. The motivation for the last 
two arguments of the scaling function in equation ([5]) comes for the behavior of the 
correlation length either with time, £ ~ t 1 ^, or with temperature, £ ~ \T — T c \~ 1 ' 1 ' . 
Letting b = t 1 ^ in equation (j3j), we obtain 

C(f, t) = t~ v/z C(r/t 1/z , \T - T c \t 1/vz ) (6) 

The algebraic prefactor corresponds to the critical behavior while the scaling function 
includes all corrections to it. The characteristic time 

T ~ C ~ \T - T c \~ uz (7) 

appears as the relaxation time of the system. In the following, we are interested only 
in autocorrelation functions, i.e. r = 0. Moreover, we expect an exponential decay of 
the scaling function Citjr) in the paramagnetic phase. Therefore, the autocorrelation 
function generally reads at equilibrium 

C(t, T) ~ e ~^JT (8) 

Two modifications need to be made to apply this hypothesis to the AFIM and the 
FFIM. First, as already mentioned, the existence of topological defects slows down the 



% The exponent in the stretched-exponential law is usually denoted f3. We will use the notation k 
instead in order to avoid confusion with the inverse temperature ft = l/fegT. 
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motion of domain walls and thus the growth of the correlation length, i.e. £ ~ (t/ hit) 1 / 2 . 
A logarithmic correction has to be included in the scaling hypothesis ([6]) 

c(o,t) = r^c(o, \t - T c \{t/^t) 1/vz ) ~ (9) 

The algebraic decay in front of the scaling function is not affected by logarithmic 
corrections because it describes the critical behavior for which topological defects are 
paired. In contradistinction, the scaling function C corresponds to the deviation to 
this behavior caused by non-vanishing scaling fields. The latter brings the system 
out-of-criticality, i.e. into the paramagnetic phase where free topological defects are 
encountered. C should thus involve logarithmic corrections. Note that the same kind 
of behavior was assumed during aging ([3]). The second modification concerns the 
relaxation time. Since the correlation length is known [3j [10] to diverge exponentially 
with temperature, i.e. £ ~ e 2 / T and not algebraically, equation (J2J) has to be replaced 
by 

T „ ~ e 2z / T . (10) 
3. Relaxation of the AFIM 

We have studied the AFIM and FFIM by means of large-scale Monte Carlo simulations 
for a two-dimensional lattice with 192 x 192 sites. The dynamics is the heat-bath local 
Markovian process introduced by Glauber [5] . We shall consider first the AFIM. Inverse 
temperatures (3 = 1/fceT in the range [0.75; 3.00] have been considered. Data have 
been averaged over 30, 000 independent histories. Error bars on the average correlation 
C(t, s) have been estimated at each times t and s as the standard deviation of the data 
produced during the different independent histories. 

To check that the system had thermalized when the measurements were started, 
we monitored the two-time spin-sin correlation functions C (t, s) which are expected to 
depend only on t — s in the stationary state. On figure [[J spin-spin autocorrelation 
functions of the AFIM are plotted at different temperatures with respect to t — s. For 
sufficiently large values of s, one can observe the collapse of the curves, indicating that 
the system has reached equilibrium. In the example of the data presented in figure [U 
one is led to the conclusion that below the inverse temperature (3 = 1/fcgT = 1.75, 
equilibrium has already been reached at time s = 1000. At (3 = 1.75, equilibrium is 
reached only at time s = 2000 while for larger values of (3, no collapse is observed 
indicating that equilibrium is not reached yet at time s = 4500. For this reason, only 
inverse temperatures (3 < 1.75 with s > 2000 will be considered in the following. 

We have tested three possible scenarii: exponential decay of the spin-spin 
autocorrelation functions (jHJ), exponential decay with logarithmic corrections ([9]), and 
stretched exponential (jlj. On figure [21 the scaling function C(t,s)(t — s) 1 / 4 is plotted 



Relaxation at finite temperature in Fully- Frustrated Ising Models 5 




fi =0 75 2000 4000 6000 8000 



Figure 1. Relaxation of the spin-spin correlation function C(t, s) of the AFIM versus 
t — s at six different inverse temperatures (3 = 0.75, 1.25, 1.50, 1.65, 1.75, and 2.00. 
The different curves correspond to different times s. Error bars are represented but 
they are hardly visible on the figure because they are of order 0(1O -5 ), i.e. much 
smaller than the symbols. 



versus t — s with a waiting time s = 2000 (ensuring equilibration as discussed above). 
In the scenario (jHJ), this function is expected to decay exponentially. This behavior is 
indeed observed over a large range of times t — s for all inverse temperatures /3 < 1.75. 
We estimated the relaxation time by interpolation of the scaling function as e~'/ r with 
a sliding interpolation window. As can be seen in the inset of figure [2} the dependence 
of t on (3 is well described by an exponential, as expected for the AFIM (see equation 
[T0l) . A fit gives r ~ e 4 31 ( 2 )/ 3 ; a behavior which is close to the expected one ( fTUl) . thought 
significantly outside error bars. 

On figure El the second scenario (J9j) is tested. The scaling function C(t, s)(t — s) 1 / 4 
is plotted versus (t — s)/ ln(t — s) . The interpolation with an exponential gives the values 
represented in the inset. As before, the relaxation time behaves exponentially with the 
inverse temperature. A fit gives the law behavior which is much closer 

to ( |T0l) than without logarithmic corrections, even though still outside error bars. 

Finally, the stretched exponential scenario is tested on figure HI The expected 
power-law behavior in the long-time regime (according to equation |4]) is observed only 
for a narrow interval of times, much narrower than in the two previous scenarii. We 
have nevertheless interpolated the data with equation (jlj) to extract the relaxation time. 
The result is plotted in the inset of figure HI An exponential growth is observed but 
with a quite different factor: r ~ e 2 - 79 ( 2 )/ 3 . 
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Figure 2. Scaling function C(t,s)(t — s) 1 / 4 at time s = 2000 versus t — s for the 
AFIM at different inverse temperatures (3 — 0.75, 1.25, 1.50, 1.65, and 1.75 (from left 
to right). In the inset, the relaxation time r obtained by interpolation of the scaling 
function as e~^~ s ^ T is plotted versus the inverse temperature (3. The straight line is 
the interpolated behavior r ~ e 4 - 31 ( 2 )/ 3 . 
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Figure 3. Scaling function C(t, s)(t - s) 1 / 4 at time s = 2000 versus (t — s)/ ln(t — s) 
for the AFIM at different inverse temperatures (3 = 0.75, 1.25, 1.50, 1.65, and 1.75 
(from from left to right). In the inset, the relaxation time r obtained by interpolation 
of the scaling function as e -(*-s)/rM*-s) i s 

plotted versus the inverse temperature (3. 
The straight line is the interpolated behavior r ~ e 3 92 ( 2 )^. 
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Figure 4. Scaling function — hxC(t, s) at time s — 2000 versus t — s for the AFIM at 
different inverse temperatures (3 = 0.75, 1.25, 1.50, 1.65 and 1.75 (from top to bottom). 
In the inset, the relaxation time r obtained by interpolation of the correlation function 

as e - ^ - ^) is plotted versus the inverse temperature f3. The straight line is the 
interpolated behavior r ~ e 2 - 79 ( 2 )P . 

4. Relaxation of the FFIM 

The FFIM has been studied in the zig-zag bond configuration for a 192 x 192 square 
lattice. We restricted ourselves to inverse temperatures smaller or equal to (3 = 1. As a 
consequence, smaller times and fewer histories were necessary. Monte Carlo data have 
been averaged over 1,000 independent histories. On figure El spin-spin autocorrelation 
functions C (t, s) are presented for different temperatures. The collapse of the different 
curves indicates that thermalization is achieved very rapidly. Like for the AFIM, we have 
compared the three scenarii: exponential decay of the spin-spin autocorrelation functions 
OH]), exponential decay with logarithmic corrections ([9]), and stretched exponential (J1J). 

First, we present the test of the stretched-exponential scenario for the FFIM 
(figure [6]). The logarithm — lnC(i, s) displays a behavior similar to the AFIM. At first 
sight, the power- law behavior that is expected according to equation (jlj) seems to be 
observed over a larger range of times t—s than for the AFIM but one should keep in mind 
that the inverse temperatures are different on figures 0] and El We have estimated the 
exponent k of the stretched exponential. Our data confirm the general trend observed in 
references [21 H]. The exponent k indeed strongly depends on the temperature. However, 
we do not recover a purely exponential decay above the percolation temperature of the 
Fortuin-Kasteleyn clusters. We do not observe two regimes (k = 1 for T > T p and k 
decreasing below T p ) as previously found but instead, an exponent slowly approaching 
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Figure 5. Relaxation of the spin-spin correlation function C(t, s) of the FFIM versus 
t — s at six different inverse temperatures f3 = 0.50, 0.60, 0.70, 0.80, 0.90 and 1.00. 
The different colors correspond to different times s as indicated in the legend. 




Figure 6. Logarithm — \nC(t, s) of the spin-spin correlation of the FFIM at different 
inverse temperatures (3 = 0.50, 0.60, 0.70, 0.80, 0.90, and 1.00 (from top to bottom). 
Only the time s = 160 is presented. In the inset, the relaxation time r obtained 

by interpolation of the correlation function as e~(~ ~) is plotted versus the inverse 
temperature /?. The two straight lines are the interpolated behaviors r ~ e 3 - 60 ^ 5 ^ for 
P > 0.6 and r - e 2 50 ^)' 3 for @ < 0.6. 
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Figure 7. Exponent k of the stretched exponential versus the inverse temperature (3 
for both the AFIM (black) and the FFIM (red). 



a value k = 1 (see figure CJ). The discrepancy may be explained by the larger lattice 
size used in this study (L = 192 instead of L = 64) and by the difficulty to identify 
a sufficiently large power-law regime, especially at low temperature. As shown in the 
inset of figure |H1 the relaxation time r does grow exponentially over the whole range 
of inverse temperatures /3 considered, in contradistinction to the theoretical prediction 
( TTUj) . However, our data are compatible with two distinct regimes of exponential growth 
with a prefactor estimated to be 3.60(5) for f3 > 0.6 and 2.50(2) for (3 < 0.6 (note that 
A ~ 0.59). 

On figure EJ the second scenario is tested. As expected, the scaling function 
C(t, s)(t — s) 1 / 4 displays an exponential decay with t — s over a large range of times t — s 
for different temperatures. However, the relaxation time does not grow exponentially 
over the whole range of temperatures considered (inset of figure [H])- Again, two regimes 
can be distinguished: the prefactor in the exponential ([TO]) is estimated to be 4.78(2) for 
f3 > 0.6 and 3.72(5) for j3 < 0.6. These values should be compared with the theoretical 
prediction 4. 

On figure El the third scenario involving logarithmic corrections is tested. The 
scaling function C(t, s)(t — s) 1//4 decay exponentially with (t — s)/ ln(t — s) over a large 
range of times t — s for different temperatures. The relaxation time r now displays an 
exponential growth over the whole range of temperatures considered and our estimate 
4.097(26) of the prefactor in the argument of the exponential is much closer to the 
theoretical prediction 4 according to equation ( ITOl) than without logarithmic corrections. 
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Figure 8. Scaling function C(t, s)(t — s) 1 / 4 at time s = 160 versus t — s for the FFIM 
at different inverse temperatures f3 = 0.50, 0.60, 0.70, 0.80, 0.90, and 1.00 (from left 
to right). In the inset, the relaxation time r obtained by interpolation of the scaling 
function as e~^~ s ^ T is plotted versus the inverse temperature fi. The two straight 
lines are the interpolated behaviors t ~ e 4 - 78 ( 2 )^ f or p > 0.6 and r ~ e 3.72(5)/3 f or 
P < 0.6. 



Model 


Stretched exponential 


Exponential 


Log. corrections 


AFIM 


2.79(2) 


4.31(2) 


3.92(2) 


FFIM 


3.60(5) {(3 > 0.6) 


4.78(2) (p > 0.6) 


4.097(26) 



Table 1. Prefactor in the argument of the exponential growth (TTU1) of the relaxation 
time for the two models AFIM and FFIM in the different studied scenarii. The 
theoretical prediction is 4. 

5. Conclusions 

We have analyzed the decay of the equilibrium two-time spin-spin correlation functions 
in the paramagnetic phase against three scenarii: exponential decay of the spin-spin 
autocorrelation functions OH]), exponential decay with logarithmic corrections (Q, and 
stretched exponential (T4|). To distinguish between them, we tested the thermal behavior 
( fTOj) of the relaxation time in the different scenarii. Our results are summarized in table 
[TJ The best agreement with the theoretical prediction 4 is obtained with an exponential 
decay with logarithmic corrections ([9]) for both the AFIM and the FFIM. Error bars 
are unfortunately smaller than the deviation from the theoretical prediction. A great 
care has been taken in the computation of error bars (from the data production to the 
analysis) so we can only invoke finite-size effects or other systematic deviations. The 
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Figure 9. Scaling function C(t, s)(t — s) 1 / 4 at time s = 1600 versus (t — s)/ \n(t — s) for 
the AFIM at different inverse temperatures /3 = 0.50, 0.60, 0.70, 0.80, 0.90, and 1.00 
(from bottom to top). In the inset, the relaxation time r obtained by interpolation 
of the scaling function as e -(*-s)/Tln(t-s) jg 

plotted versus the inverse temperature j3. 
The straight line is the interpolated behavior r ~ e 4 ' 097 ( 26 )^. 

deviation for the two other scenarii being much larger, we do not expect them to get 
closer to 4 than the exponential decay with logarithmic corrections. These results bring 
further evidences of the existence of topological defects in the paramagnetic phase of 
the AFIM and the FFIM. 
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